arXiv:1508.00630vl [math.NA] 4 Aug 2015 
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Abstract 

We develop and analyze strategies to couple the discontinuous Petrov-Galerkin method 
with optimal test functions to (i) least-squares boundary elements and (ii) various variants 
of standard Galerkin boundary elements. Essential feature of our methods is that, despite 
the use of boundary integral equations, optimal test functions have to be computed only 
locally. We apply our findings to a standard transmission problem in full space and present 
numerical experiments to validate our theory. 
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1 Introduction 

In its current form, the discontinuous Petrov-Galerkin (DPG) method with optimal test func¬ 
tions has been introduced by Demkowicz and Gopalakrishnan mm- Principal objective is to 
guarantee uniform stability of the discrete scheme by the use of specific optimal test functions, 
in particular for singularly perturbed mum and wave problems GS3H3I22]. 

Essential feature of this method is that optimal test functions be calculated (approximatively) 
locally on elements. This is well understood for partial differential equations on bounded domains. 
When dealing with problems in unbounded domains (like transmission problems) boundary inte¬ 
gral equations are a natural way to deal with the exterior part. Their numerical analysis usually 
requires the use of trace spaces, which are Sobolev spaces of non-integer orders. As the integral 
operators themselves, corresponding norms are non-local. That is, they cannot be represented 
equivalently as broken norms over elements and, therefore, there are no straightforward efficient 
techniques for the calculation of optimal test functions in these cases. Perhaps surprisingly, in 
[26i i 2Tj we found an ultra-weak variational formulation for hypersingular operators that is well 
posed in integer-order Sobolev spaces as long as the underlying polygon (in two dimensions) or 
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surface (in three dimensions) is closed. Then, also the DPG framework with optimal test func¬ 
tions goes through without complications (this is different for open curves and surfaces). Still, 
the underlying formulation involves boundary integral operators. This implies that optimal test 
functions cannot be calculated locally, and this seems natural for global operators. 

Now, returning to transmission problems with variational formulations comprising boundary 
integral operators, application of DPG technology faces the problem of how to efficiently calculate 
optimal test functions. In [[25] we propose to globally apply this technique without specific 
separation between differential and integral operators. It turns out that, as expected, the integral 
operators generate global effects for optimal test functions, but only at a strip along the interface 
(on which the integral operators live). In this way, one partially looses locality in the calculation 
of optimal test functions. 

In this paper, we pursue the idea of splitting approximations for partial differential operators 
from those of boundary integral operators. This is very natural when coupling finite elements 
with the boundary element method (BEM). Indeed, there is a long list of references and we only 
cite [28i nni eta go to name a few classical ones, and refer to more recent developments [MIESIEI] 
for detailed discussions. However, when trying to use similar coupling techniques of DPG with 
BEM, there is a certain conflict of frameworks. DPG-analysis is based upon arguments from 
functional analysis (even at the discrete level) and the assignment of optimal test functions to 
any approximating basis function. This generates square systems and complicates the coupling 
with Galerkin boundary elements where some of the approximating finite element basis functions 
are tested with traditional basis functions from other discrete test spaces. 

Noting that the DPG method with optimal test functions is in fact a least-squares method, 
it seems most natural to couple it with a least-squares boundary element method. Previously, 
least-squares boundary elements have been coupled with finite elements in [20 s 3ly EQ], and also in 
combination with the DPG method they perform well, as we will show. Additionally, we propose 
three coupling methods with Galerkin boundary elements, one with hypersingular operator, one 
with weakly singular operator, and another one with a combination of both. 

Differently from traditional coupling methods, where some of the Cauchy data are represented 
directly via boundary integral operators, our boundary element equations work as constraints 
for the DPG scheme. In this way, ellipticity of the resulting bilinear form can be achieved and 
the well-posedness of the discretization follows from that of the continuous formulation. We 
note, however, that in the case of the Galerkin variants, proofs of ellipticity require weighting 
of the DPG part with a positive number whose size is in principle unknown. In our numerical 
experiments we always choose the constant one and observe perfect results. Numerical tests for 
our model problems (and the specific domain under consideration) indicate that only a small 
weighting of the order of 1/3 or less results in a stiffness matrix with symmetric part which is 
not positive definite. 

Essential advantage of all our coupling methods is that optimal test functions have to be cal¬ 
culated only for terms involving differential operators. In this way, locality of the corresponding 
trial-to-test operator (that maps ansatz functions to optimal test functions) is fully maintained 
and DPG technology is applicable without interference of non-local operators. In fact, we expect 
that in this way DPG formulations, that are specifically designed for singularly perturbed and 
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wave problems, can be coupled to boundary elements. On the other hand, the analysis of two 
of our methods (least-squares BEM and hypersingular operator) is based upon the fact that 
the kernel of the interior formulation is strongly related to solutions of the exterior part of the 
problem. In fact, Calderon boundary integral operators (used for our coupling schemes) map 
kernel functions of the interior problem to components of their Cauchy data. This relation is 
maintained whenever the interior and exterior differential operators of the transmission problem 
are identical. At the moment it is unclear whether our analysis of those two methods can be 
extended to the case of different interior and exterior problems. 

To prove ellipticity of the coupling variant with weakly singular operator and Galerkin ap¬ 
proximation, we apply a recent technique from Sayas and Steinbach |341 135] . Here, only gentle 
relations between Cauchy data of interior and exterior problems are needed, and the proof ex¬ 
tends to interior operators of the form —divA(x)Vrr with minimum eigenvalue of A strictly 
bounded below by 1/4. 

The analysis of our fourth coupling scheme with a combination of hypersingular and weakly 
singular operators, also does not rely on a specific relation of Cauchy data. It is the coupling 
variant which is relatively straightforwardly applicable to more general problems combining, e.g., 
singularly perturbed PDEs on bounded domains with linear, homogeneous PDEs and constant 
coefficients in the exterior. Contrary to our analysis of the coupling with weakly singular oper¬ 
ator, such an extension would not need a condition on the minimum eigenvalue of the interior 
differential operator. 

The remainder of this paper is as follows. In Section [2] we present the model problem, recall 
some Sobolev spaces and norms (Section 12.2p . revisit properties of boundary integral operators 
(Section 12.31) . present our three coupled schemes (Section I2.5p . and state their well-posedness 
and quasi-optimal convergence (Theorems 01010) . The presentation of the schemes is first done 
in an abstract way (in Section m- We then recall an ultra-weak variational formulation of the 
interior part of the problem (Section 12. 4h and deal with the three coupling schemes, respectively, 
in Subsections 12.5.11 12.5.21 and 12.5.31 Further technical details and proofs of the main theorems 
are given in Section 0J In Section 0] we present several numerical examples. We also give details 
on the implementation of our least-squares boundary elements (Section 14. ID and present a simple 
a posteriori error estimator (Section 14.2(1 that serves as a common error bound for comparison. 

2 Mathematical setting and main results 

Let C M d , cL £ {2,3}, be a bounded, simply connected Lipschitz domain with boundary T, 
fi c := M. d \ fi, and normal vector on T pointing in direction of H c . We consider the following 
model transmission problem: given / £ ^(tl), uq £ iL 1 / 2 (T), </o £ T), find u £ iL 1 (I2) 
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and u c € H 2 oc (Q c ) such that 


—A u = f in fl, 

A u c = 0 in I2 C , 
u — u c = uq on T, 
d 

g^( u - u ) = fa on r, 

u c (x ) = 0 ( |x| _1 ) as \x 


(la) 

(lb) 

(lc) 

(l d) 

oo. (le) 


Here, H^ oc (fl c ) denote standard Sobolev spaces with trace space H l ^ 2 (T) and its dual 

h -v 2 (r) , see Section [2721 for precise definitions. For d = 2, we assume in addition that f (} f + 
J r 00 = 0 (which ensures the correct behavior ([Tell at infinity). 


2.1 Abstract coupling framework 

Our aim is to use a DPG formulation for the interior part of problem ((TJ), and to couple it with 
boundary integral equation(s) to take the exterior part into account. We will propose four coupled 
schemes whose settings are formally identical. In this section we present the corresponding 
abstract framework. 

Let U and V be, respectively, a reflexive Banach space with norm || • ||u and a Hilbert space 
with inner product (•, -)y and norm || • ||y. For a given bilinear form &(•, •) : U x V —> M and 
linear functional Ly G V' we consider the variational formulation 

u G U : b(u,v ) = Ly(v) Mv G V. (2) 

We will also make use of the operator B : U —> V' induced by the bilinear form b. In our case, 
problem (J2]) itself has no unique solution and will be coupled with a variational relation 

u £U : c(u,w) = Ljj(w ) Vic € U. (3) 

Here, c(-, •) : U x U —> R is another bilinear form and Ljj € U' a linear functional. Essential 
feature in the discretization of m is the use of so-called optimal test functions. They are 
generated through the weighted trial-to-test operator Qp : U —>■ V defined by 

(@pu,v)y = /3b(u,v) \/v£V (/3 > 0). (4) 

Using this operator for a constant /3 > 0, we combine the relations (J2J) and ([3]) to the following 
variational formulation: 

u£U: b(u,Qpw) + c(u : w) = Ly(Qpw) + Lu(w) \/w £ U. (5) 

Note that in this DPG-type formulation, the trial-to-test operator makes use of the bilinear form 
&(•, •) exclusively. In our specific schemes introduced below, this allows for maintaining locality 
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of the application of 0^ while the bilinear form c(-, •) comprises non-local (boundary integral) 
operators. 

For a finite-dimensional approximation space Uh p C U, the discrete scheme then is 

u hp <E U hp : b(u hp , Qpw) + c{u hp , w ) = L v (®pw) + L v (w) Vw € U hp . (6) 

The following lemma states well-posedness of ([5]) and best-approximation property of its discrete 
form ([6]). It amounts to using the Lax-Milgram lemma (in the case of Hilbert spaces) or Babuska- 
Brezzi theory (in the case of Banach spaces) and Cea’s lemma for elliptic variational formulations. 
Of course, this is a standard procedure in least-squares analysis, see, e.g., [4j. 

Lemma 1. Suppose that b : U x V —> M and c : U x U —>• M are bounded bilinear forms with 
bounds Cb and C c , respectively. Furthermore, assume that, for given fd > 0, there is a > 0 such 
that 


a\\w\\u < b(w,Q/ 3 w) + c(w,w) \/w £ U, 

that is, the bilinear form on the left-hand side of (JSJ) is U-elliptic. Then, there is a unique 
solution u of (f5|) with 


u 


\u < —(PCb\\Lv\\v' + \\Lu\\u')- 

a 

In addition, there is a unique solution u^p of (|6|), and it satisfies 

pc 2 b + C, 


W - u hp \\u < 


a 


inf ||m — 

wGU hp 


Proof. The statements are immediate consequences of the Lax-Milgram lemma (or Babuska- 
Brezzi theory) and Cea’s lemma. In fact, ellipticity of the bilinear form holds by assumption 
and boundedness of the bilinear form and linear functionals follows by bounding ||© ( gt«||y < 
/3C&||iu||[/ for all w € U. To this end, note that 0^ = /3Jy^B with isometric Riesz map 
Jy : V —> V'. To show the best approximation property we also make use of the orthogonality 
b{u - u hp , Qpw) + c(u - u hpi w) = 0 Vru <E U hp . □ 

In Section 12.41 we will specify ([2]) as an ultra-weak variational formulation of the interior part 
of m i an d in Sections 12.5.11 12.5.21 12.5.31 we will present three possibilities of selecting ([3]) as 
boundary integral equations. Before doing so we introduce some Sobolev spaces and norms in 
Section 12.21 and recall properties of boundary integral operators in Section 12.31 


2.2 Sobolev spaces and norms 

For a Lipschitz domain u C W l we use the standard Sobolev spaces L 2 (w), H l { oj), Hq(oj) (the 
index 0 denotes vanishing trace on the boundary of w), H(div,cj), Ho(div,w). Vector-valued 
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spaces and functions will be denoted by bold symbols. Denoting by 7 ^ the trace operator acting 
on i7 1 (w), we define trace spaces 

H l / 2 (duj) := { 70 ju; u G H 1 ( w)} and its dual H~ l t 2 (du) := (iL 1 // 2 (<9u;)) / 

and use the canonical norms. Here, duality is understood with respect to L^dui) as a pivot 
space, i.e., using the extended L 2 (du) inner product (-,-) 9uj . The inner product will be 

denoted by Let T denote a disjoint partition of H into open Lipschitz sets T G T, i.e., 

U^eTT = H. The set of all boundaries of all elements forms the skeleton S := {dT | T G T}. 
By n m we mean the outer normal vector on dM for a Lipschitz set M. On a partition T we use 
product spaces H ] (T) and H(div,T) and equip them with corresponding product norms. The 
symbols V 7 - and div^ denote, respectively, the T-piecewise gradient and divergence operators. 
On the skeleton S of T we introduce the trace spaces 

H 1/2 (S) := ju G U Te rH 1/2 (dT)- 3w G H 1 ^) such that u\ 9T = w\ dT VT G t}, 

H~ l / 2 (S) := jcf G HT&TH~ l / 2 {dT)] 3q G H(div, H) such that Sj dt = (q ' n T) 1 8T VT G t|. 

These spaces are equipped with the norms 

II“IIj?i/ 2 ( 5 ) : = inf {||wHj/i(o); w G iL 1 (H) such that u\ 9T = w\ dT VT G T) , (7a) 

11511^-1/2(5) := inf {||<?|| H (div,0); Q e H(div , fi) such that a\ 9T = (q ■ n T )\ 9T VT G T} . (7b) 

For our analysis we will also need the space 

:= {ueH^iSy, u|r = 0}. 

For functions u G H 1 ^ 2 (S), a G H~ l / 2 (S) (which are elements of product spaces with components 
u\ 9 T G iL 1 - / 2 ((?T) and d\ 9 t G H~ 1 / 2 {dT), respectively, for T G T), and r G H(div,T), v G 
H l (T) we use the notation 

(u,r • n ) 5 := ^(u\ 9T ,r ■ n T ) 9T , (a,v) s ■= ^ (°\dT,v)dT- 
TeT Ter 

Furthermore, throughout the paper, suprema are taken over sets excluding the null element, and 
the notation A < B is used to say that A < C ■ B with a constant C > 0 which does not depend 
on any quantities of interest. Correspondingly, the notation A > B is used, and A ~ B means 
that A < B and B < A. 


2.3 Boundary integral operators 


The exterior part of problem (JT]) will be dealt with by boundary integral operators. To this end 
we need some further definitions. The fundamental solution of the Laplacian is 


G(z) 


_ ^7T l0 & \ Z 

J_ 1 

4 tt HI 
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(5 = 2), 
(5 = 3), 


and the corresponding single layer and double layer potentials are 

V(f>{x) := J G(x — y)(j>(y) ds y , Kv(x) := j d n ^ y )G(x - y)v(y) ds y , x <5 M d \ T. 

Application of the trace operator jq and the normal derivative <9 nn gives rise to the three 
boundary integral operators 

V:= 7 cV, /C:=1/2 + 7q /C, W := -8 nn K. (8) 

They are the single layer, double layer, and hypersingular operators, respectively. The adjoint 
operator of 1C is denoted by 1C'. These operators are linear and bounded as mappings V : 

-7 tf 1/2 (T), K : tf 1/2 (T) -7 H 1 ' 2 ^), K! : iU 1/2 (r) H^/ 2 (T), and W : H l / 2 {V) ->• 

i4 _1 / 2 (r). We note that here holds ker(W) = span{l}. For (u,4>) € Ff 1//2 (T) x H~ l / 2 (T) let us 
define the operators 

¥(u, </>) := Vcj) + (1/2 - K)u, W {u, </>) := Wu + (1/2 + 

By the boundedness of V, 1C, 1C', and W, it follows that 

¥ : H 1/2 (T) x H~ 1/2 (T) ->■ H 1/2 (T), 

W : H 1/2 (T) x H~ 1/2 (T ) -> FT 1/2 (r) 

are bounded. Furthermore, for the exterior Cauchy data of the harmonic function u c £ H^ oc (Cl c ) 
(see fllbli ) we have the so-called Calderon system 

¥(u c | r ,5 nn u c ) =0, (10a) 

W(u c | r , d nn u c ) = 0. (10b) 

For details and proofs we refer to classical references, e.g. [9U32U27]. 

2.4 Ultra-weak finite element part 

In this section we recall an ultra-weak variational formulation of the interior part of the trans¬ 
mission problem ©■ This is taken from [13] and corresponds to problem ([2]) of the abstract 
framework in Section 12.11 

To this end, let T be a partition of Cl with skeleton S. Corresponding to T we define 
fractional-order spaces H 1 ^ 2 (<S) and H~ 1 ^ 2 {S) as in Section T2.21 and select 

U := L 2 (tt) x L 2 (fi) x iU/ 2 (S) x H~^ 2 (S), V := H\T) x H(div,T). 

The broken test space V is provided with the canonical inner product denoted (as previously) 
by (‘) ')v- Furthermore, we define the bilinear and linear forms b and Ly by 

b(u,v) := (u, div r r)n + (er, V r v + t)q - (u,t ■ n) s - (d,v) s , L v (v) := (f ,v)n (11) 
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for u = (u, cr , u, a) £ U, v = ( v , t) £ V. 

With this setting, © is our ultra-weak variational formulation of the interior part dial) of 
the transmission problem. For reference, we explicitly specify its strong form 

u := (u,cr,u, a) £ U : Bu = Ly. (12) 

The weighted trial-to-test operator 0^ is exactly as in Q. 

Contrary to [T5] (where the model problem considers a homogeneous Dirichlet boundary 
condition), our transmission problem uses a space U whose third component H 1 ^ 2 (S) does not 
incorporate such a condition on T. We therefore anticipate that the operator B : U —> V' has a 
non-trivial kernel whereas it is known that 

B: U 0 := L 2 (Q) x L 2 (D) x H^ 2 (S) X H~ 1/2 (S ) -> V' 

is invertible: 

Lemma 2. The bilinear form b : U X V —> M is bounded with a constant independent of T■ 
Furthermore, there holds 

sup V ^ > ||mo||i/ f or a U u o € U 0 , and 

v&v ll^llv 

sup ^, 1 1X0 1 ^ > 0 for all v £V\ {0}. 
uoeUo \\ u o\\u 

Proof. The boundedness of b is immediate by the definition of the norms. The remaining results 
follow from [13]. More specifically, by m Theorem 4.2] there holds the converse inf-sup condition 
II-BuIIj// > HvHy Vd £ V, and b is non-degenerate in the [/o-component by [13l Lemma 4.1]. The 
statement then follows by the Babuska-Brezzi theory. □ 

Obviously, the non-uniqueness of (11211 is due to missing boundary conditions. They will be 
incorporated by Cauchy-data relations that stem from the transmission conditions ([Tc]) . dldl) and 
the exterior problem dlbl) . (flcl) . For the handling of these data we need the following restriction 
operators. 

7 : U — > # 1 / 2 (r) x Ff - 1 / 2 (r), 7 (u,cr,u,a) := (u|r,?|r) (Cauchy data), 

70 : U H ^(r), 7 o (u,<r,u,a) := u|r (Dirichlet trace), 

7 n : U —> H~ l / 2 (T), 7 n (ri, cr, u, a) := u|r (Neumann trace). 

The boundedness of these operators is immediate. 

Lemma 3. The operators 7 , 7 n and 70 are bounded. 

Proof. By definition of Ff 1 // 2 (T) and its norm, ||M||^i/ 2 ( r ) < ||w||j/i(n) for any u £ iL 1 (D) with 
u = u on r. That is, ||u||^i/ 2 / r ) < ||u||^i/ 2 ( 5 ) by the definition of the latter norm. Analogously, 
making use of the boundedness of the normal component trace from H(div , fl) to iW 1 // 2 (T), we 
obtain 





2.5 Four coupling methods 

In Section 12.41 we have fixed the interior part of the variational formulation for the transmis¬ 
sion problem ([T]). This corresponds to the abstract form Q. It remains to add an equation 
([3]) to incorporate the exterior part of the problem and transmission conditions. In the next 
subsection we consider a boundary integral equation of least-squares type for Q, and in Sub¬ 
sections 12.5.2112.5.3112.5.41 we study three Galerkin variants, with hypersingular operator, single 
layer operator, and both operators, respectively (all with additional rank-one terms). 

2.5.1 Coupling with least-squares boundary elements 

We couple equation (1121) with the exterior problem by using the boundary integral equation (llOal) 
and jump relations m, d- With the operator V this yields the relation 

V(7«) = V(u 0 , 4>o) in i7 1/2 (T). (13) 

Denoting the inner product in H l / 2 (T) by (•, •) H i/ 2 ^y a least-squares formulation of (1131) is 

(V(7w),V(7w)) ff i /2(r) = (V(n o ,0o),V(7m)) H i/ 2(r) Vrn € U. 

According to the abstract framework from Section 12.11 we define 

c(u,w) :=c LS {u,w ) := (Y(-fu) ,Y(-/w)} H i /2{r) , 

Lu{w) := L ls (w) := (Y(u 0 ,(po) ,Y('yw)) Hl / 2 ^ r y 

Then our DPG formulation coupled with a least-squares boundary integral equation is (f5|) : 

u£U: b(u,@i 3 w) + cls(u,w) = LviQpw) + Lls(w) Vrn €U. (14) 

The bilinear form b and functional Ly have been defined in Section 12.41 The coupled scheme 
(1141) is equivalent to the least-squares formulation 

u = (u,a,u,a) G U : u = argmin Jp{w\ /, uq, </>o) (15) 

wGU 


with 

:= P\\Bu-L v \\y, + ||V(7«) - Y(u 0 , </>o)||^i/ 2 (r) - 

The corresponding discretization is ([HI) with U replaced by a finite-dimensional subspace U^p, 
cf. (J6]) , or equivalently 


u hp € U hp : u hp = argmin Jp(w;f,u 0 , )■ (16) 

Principal advantage of this scheme over the following ones is that, for scaling parameter j3 = 1, 
the resulting bilinear form is known to be elliptic so that the formulation is well posed and its 
discretization converges quasi-optimally. 
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Theorem 4. Let (3 = 1. Problems © and (usd are equivalent and uniquely solvable. More 
precisely, let (u, u c ) be the solution to © and define u := (u,cr,u,i t) with a := Vu in Ll, 
u := u\s (element-wise), and a := er • ny on dT for any T G T. Then, u satisfies (j!5|) . On 
the other hand, if u = (u,a,u,a) solves (fT5]l . then (u,u c ) with u c := /C(u|r — uq) — V(<r|r — 4>o) 
solves ©• 

Furthermore, for a finite-dimensional subspace Uh p C U, there is a unique solution Uh p G Uh p 
of (USD , and there holds the best approximation property 


\\u-u hp \\u< inf \\u-w\\u 
w£U hp 

with hidden constant independent of T and the polynomial degrees in Uh p . 
A proof of this theorem will be given in Section 13.11 


Remark 5. (i) In contrast to the method proposed in J25]/ , the trial-to-test operator 0 := ©i used 
here can be implemented without considering any boundary integral operator. Therefore, having 
defined the space V as a product space of broken Sobolev spaces of integer orders, the action of 
0 is completely local and can be approximated in a standard way known from DPG methods, 
cf. 

(ii) The only difference of our trial-to-test operator with that of standard DPG schemes is 
that, in our case, the bilinear form on the right-hand side defining 0 has a kernel. This does not 
affect the implementation. Indeed, if there is w G Uh p n ker B then Qw = 0 and equation ([HI) 
reduces to 


(Y ( r yuh p ) , r yo w ) H i/2(r) = (V(u o ,0o),7o^)ifi/2(r), 

cf. Lemmas \Tl\ and 1 12[ 

(Hi) The boundary element least-squares term (V(7-), V( 7 -))#i/ 2 (p) is of a standard type and 
has to be approximated by using an appropriate preconditioner instead of the H 1 ^ 2 (r) -inner prod- 

; 30) for similar techniques. 


uct. We give precise details in Section f.l and refer to f20[ 131 

(iv) It is also possible to define a least-squares coupling with equation (UObl) instead of (1 lOajl . 
Due the fact that W(1,0) = 0 , a (rank-one) stabilization term has to be added to ensure unique¬ 
ness. Here, the same term as for the coupling with Galerkin boundary elements can be used (see 
Subsections \2.5.21 \2.5.3\ below). 


2.5.2 Coupling with Galerkin boundary elements (hypersingular operator) 

We proceed as in Subsection 12.5.11 but couple equation (11211 with a variational form of (11 Obi) 
including transmission conditions instead of a least-squares form of (HOa.p . that is, 

(W(yu) ,7 0 ie)r = (W(u 0 , fio) , 7o™)r Vu; G U. 

Since W(1,0) = 0, a rank-one term has to be added for uniqueness. We therefore incorporate 
relation (HOaP combined with the transmission conditions, and define the bilinear and linear forms 

c(u,w) := c hy (u,w) := (W( 7 tt), 7 0 ir)r + (1, V( 7 u)) r (l, Y(jw)) r , 

Lu(w) := L hy (w) := (W(u 0 , </> 0 ), 7o™)r + (1 , V(u 0 , 0o))r(l , V(7it»)) r . 
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Then, instead of (fl4l) we obtain the following coupling variant of (j5|): 

u£U: b(u, Qpw) + Chy(u, w) = Ly(Qpw) + L^ y (w) \/w £ U. (17) 

Again, the bilinear form b and functional Ly are the ones from Section 12.41 Correspondingly, 
the discrete scheme is 

'U'hp £ Uhp . b(uh p , Qfjw) + Chy('U/ l p, w) — Ly(®pw) T Z/j 1 y('in) Via £ Uh p - (18) 

It turns out that the bilinear form defined by the left-hand side of (11711 is IT-elliptic only for 
sufficiently large (3. This is the reason for having introduced this parameter. The following 
theorem is our second main result. 

Theorem 6. There exists a constant (3 q > 0 which depends only on fl such that for all /3 > (3 q 
problems © and Gzi) are equivalent and uniquely solvable. More precisely, let (it, u c ) be the 
solution to © and define u := ( u,cr,u,a ) with a := Vit in u := it|,s (element-wise), and 
a := (7 ■ n t on dT for any T € T. Then, u satisfies (1171) . On the other hand, if u = (it, er, u, a) 
solves (1171) . then ( u,u c ) with u c := JC(u |r — ito) — V(a|r — <^o) solves ©. 

Furthermore, for a finite-dimensional subspace Uhp C U, there is a unique solution Uh p £ Uhp 
of (IlSp . and there holds the best approximation property 

\\u - u hp \\u < (3 inf ||it-iu||t/ 

tiJ ££/hp 

with hidden constant independent of F, Uhp, and fi (under the restriction j3 > fio). 

A proof of this theorem is given in Section 13.21 

Remark 7. (i) It can be shown that, for arbitrary f5 > 0 but with upper bound for the diameter 
of Q, the bilinear form b(-,Qp-) + Ch y (-, •) from (1171) (with slightly different stabilization term) 
satisfies both inf-sup conditions at the continuous level. However, a proof of the discrete inf-sup 
condition is unknown. We therefore stick to an analysis based on ellipticity. 

(ii) The corresponding versions of Remarks {5\ (i), (ii) apply in this case as well. 

2.5.3 Coupling with Galerkin boundary elements (single layer operator) 

We proceed again in the fashion of Subsection 12.5.11 This time we couple equation (1121) with a 
variational form of (110a[) including the transmission conditions, that is, 

(7 n w,V(7ti))r = (7n^,V(it 0 ,</ , o))r Via £ U. 

We add the same rank-one term as before to ensure uniqueness, and define the bilinear and linear 
forms 


c(u,w) := Csi(it, w) := ('j n w , V(7it)) r + (1 ,¥( 7 it)) r (l, V(7in))r, 
Lu(w) := L sl (w ) := (j n w , Y(u 0 , </> 0 ))r + (1 , V(ix 0 , </> 0 ))r(l, V(7«j)) r . 
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Then, instead of the Galerkin boundary element coupling with hypersingular integral operator 
(11711 we obtain the following coupling variant of (|5|): 

u £ U : b(u,@pw) + c s i(u,w) = Ly(Opw) + L s \(w) Mw £ U. (19) 

As before, the bilinear form b and functional Ly are the ones from Section 12.41 The discrete 
scheme is 

u-hp G Uhp . b^Ufrpj Qpw') + Csi(uh p , w) Ly(®pw) + L s i(w) \/w £ Uhp. (20) 

For this coupling variant we require that the single layer operator V is if _1 / 2 (r)-elliptic. This 
is always true in three dimensions and, in two dimensions, holds if the diameter of the domain 
is sufficiently small. Main advantage of this coupling variant is that its analysis extends to 
operators of the form —divA(x)Vu in (Hall whose minimum eigenvalues are strictly bounded 
from below by 1/4, cf. [34j [35]. (To be precise, the analysis given in Section [3.31 can be applied 
to such an operator. An improvement of the technical estimate (l29j) there, leads to the bound 
for the minimum eigenvalue.) 

Theorem 8. In two dimensions (d = 2) assume that dianr(fl) < 1. Then, the statement of 
Theorem 0 holds correspondingly for the coupling CHD and its discrete variant m- 

2.5.4 Coupling with Galerkin boundary elements (Calderon system) 

We now use both Calderon equations (1 1GH to couple the exterior part of the transmission problem 
to the interior ultra-weak formulation (this is (1121) in strong form). The equations are 

(W( 7 u ), 7oin)r = (W(u 0 , <f>o ), 7o™)r Vw £ U, 

( 7 n m,V( 7 u)) r = ( 7 ^ ,¥(u 0 , (f> 0 )) r Vm £ U. 

We add the same rank-one term as before to ensure uniqueness, and define the bilinear and linear 
forms 

c(u,w) := c ca .(u,w) := (W( 7 u), 7 0 m) r + ( 7 n m,V (7 u)) T + (1 ,V( 7 u)) r (l ,V('jw)) r , 

Lu(w) := L ca (w) := (W(u 0 , (f>o ), 7o^)r + ( 7 n™, ¥(u 0 , </ 0 ))r + (1, V(u 0 , </> 0 ))r(l, V( 7 io))r- 

Then, we obtain the following coupling variant of (0: 

u £ U : b(u,@i 3 w) + c ca (u,w) = Ly(Opw) + L ca (w) Vw £ U, (21) 

and discrete scheme 

u hp G U hp : b(u hp , Qpw) + c ca (u hp , w) = Ly(Qpw) + L ca (w ) Vm £ U hp . (22) 

Again, the bilinear form b and functional Ly are the ones from Section 12.41 For this coupling 
variant, as for the case with weakly singular operator, we need 77 _1 / 2 (r)-ellipticity of the single 
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layer operator V. As previously noted, in three dimensions this is always true whereas in two 
dimensions the domain has to be appropriately scaled. Main advantage of this coupling variant 
is that its analysis does not make use of special relations between Cauchy data of the interior 
and exterior problems (specifically, Lemma [T2l from Section [3] below is not needed). Its proof is 
therefore easier to extend to problems where the differential operator in the interior is different 
from the one in the exterior. Also, we do not require a condition on the minimum eigenvalue of 
operators of the type —div A(x)Vu for the interior problem (llajl . which is needed when extending 
the variant from Section 12.5.31 to such operators. 

Theorem 9. In two dimensions (d = 2) assume that dianr(P) < 1. Then, the statement of The¬ 
orem 0 holds correspondingly for the variational formulation m and its discrete variant (l22l) . 


3 Analysis and proofs of the main results 


In this section we provide some technical details and prove Theorems [Tj [ 6 ] and [ 8 j 
For given ip £ H 1 // 2 (r) let us define its harmonic extension 


£(<p) := u = (u, cr, u, a) 


u £ IL 1 (0) : u = ip on T, (V'u, Vtc)o = 0 Vw £ Hq (O), 
a = Vu, u = u on S, a = cr ■ n-p on dT , VT £ T. 


We can state the following properties of £. 

Lemma 10. The operator £ : H l ^ 2 {T) —> U is linear, bounded, and a right-inverse 0 / 70 . 

Proof. By construction we have 70 = <p for any <p £ H 1 / 2 ( r), showing that £ is a right- 

inverse. The continuity of £ follows from the well-posedness of the Dirichlet problem, i.e., 
IMIiUfO) 1 $ the definition of the other components of £ ip and the involved norms, 

cf. dZJ) and Lemma El □ 

Lemma 11. There holds ker B = £H l / 2 (T). 

Proof. We note that £ H l / 2 (T) C ker B follows from integration by parts. Therefore, it remains 
to show that ker B C £ H l / 2 (T). Let u £ U be arbitrary but fixed and define Uq := u — £ r )Q u - 
Then 70^0 = 0 by Lemma [TUI he., Uq £ Uq, and due to Lemma [ 2 ] it holds 

\\Bu\\ V t = ||.Btio||v" = SU P b i^i V ^ ^ II“olit/- (23) 

vev ||w||v 

Suppose that u £ ker B. From (l23l) it follows that u = £ 70 u £ £ H 1 / 2 ( r). □ 

Lemma 12. There holds ¥(7 u) = 70 u and W(ju) = 7 n u for all u £ ker B. 

Proof. Let u = (u, cr ,u,a) £ ker B be given. By Lemma [Tl] and the definition of £ we have 

^ (fjfl ^ 

A u = 0 in P, u = u|r on T, and —— = ulp on T. 

OUQ 
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Therefore, the representation formula (see |32j ) yields u = Va — K,u in 9, and application of the 
trace operator proves that u/2 = Vix — ICu on T, cf. (| 8 ]h that is, 

V (7 u) = Va + (1/2 — /C)rt|r = rt|r = 7o u - 

The relation involving W follows by a similar argument. □ 

The following norm equivalence follows by a standard compactness argument. Similar results 
for classical coupling methods can be found in [2 [T9]. 

Lemma 13. In the case d = 2 assume that diam(9) < 1. There holds 

IMItfi/ 2 (r) + IMItf—i/ 2 (r) - (yVu,u ) r + {(j>,V4>) r + |(l,V(u») r | 2 
for all ( u,4 i) G X H 1 // 2 (T), and the involved constants only depend on T. 

3.1 Proof of Theorem [4] 

It is well known that m is uniquely solvable, see, e.g., m Let (u, u c ) solve m and let u = 
(u,cr,u,a) be defined as in the theorem. Then, by construction, J\{u] /, uq, <f>o) = 0, i.e., u 
solves (1151) . It remains to show that any solution u to (I15[) is unique, and that u^p satisfies the 
quasi-optimal error estimate. Both follow from Lemma [T] 

To this end, we note that the right-hand side functionals Ly(Q l •), Lls : U —> V' are 
bounded by the boundedness of 0i (see the proof of Lemma|I|, Ly (since / G L 2 ( 9 )), V (by Q), 
and 7 (by Lemma@. It remains to show the [7-ellipticity of the bilinear form &(•, ©r) + cls(-, ■). 
Let u = (u,tT,u,a) G U be given. By Lemmas [2l 1101 andllllwe know that 

\\Bu\\yt>\\u-£^u\\u. (24) 

Furthermore, the continuity of the extension operator (see Lemma ITOl) . Lemmas [Til E21 and the 
triangle inequality show that 

\\£lou\\u< ||7o«||tfi/2 (r) = ||V( 7 £ 7 0 u)|| H i/ 2 (r) 

< l|V(7«)|| H i/ 2 ( r) + ||V(7(u-f 7ou))|| Jf i /a(r) . (25) 

Now, using again the continuity of V and 7 , estimate f|24[) yields 

II V(y(w - f7 0 w))ll^i/z^) < \\u-£” 1 qu\\tj < \\Bu\\ vl . (26) 

Finally, combining ()24[) (l26l) and the triangle inequality proves that 

\\ u \\u < II V( 7 tr)|| H i/ 2 (r) + \\Bu\\ V t yueU. (27) 

Recalling the definition of cls(v) an< ^ the relation = b(u,@iu), this proves the U- 

ellipticity of the bilinear form. □ 
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3.2 Proof of Theorem [6] 

We only have to slightly vary the proof of Theorem 3] from Section 13.11 

Integration by parts and the Calderon system (HOD show that the solution u of (JT]) also solves 
problem (|17D . It remains to prove that (117[) is uniquely solvable and that discrete approximations 
fulfill the best approximation property. Again we use Lemma [lj 

In the current case, boundedness of c(-, •) = Ch y (-, ■) follows from the mapping properties of 
the integral operators V and W (see Section 12.3D , and the boundedness of 7 and 70 by Lemma [3j 
It remains to prove t/-ellipticity of the bilinear form £>(•, 0^ •) + c(-, •). 

We use the relations from Lemma [12] and the norm equivalence ||v||^i^ — ll^ 7 t 'lli 2 (o) + 
|(1 ,u)r | 2 for v € Furthermore, let u denote the first component of £ 70 u. Integration by 

parts and the harmonicity of u in P show that (Vu, Vu)n = (7n2 ) 7o' li )r- We obtain 

c(£7 0 m,£7o«) = (yVu + (5 + K')^ n £^u , u) T + |(1 , V(7£7 0 u))r | 2 

= ( 7 n£ 7 ow,S) r + |(l,S) r | 2 = ||Vu||| 2(n) + |(l,S)r | 2 

- ll^7oM||fri ( n) ^ l|7ow||^i /2(r) > \\£~/ 0 u\\u. (28) 

By Lemmas [2] HU] and[TT]we bound 

\\ u \\u ~ ll u ~ £7o' u llr/ + II £■ lo u \\u 1$ ll^ ^l ll^" , + II £ r yo u \\u- 

Combining the last two estimates, and using the boundedness of c(-, ■) and Young’s inequality, 
we conclude that, for 5 > 0, 

\\u\\u ^ \\ Bu \\v' + l|£7o «||u < \\Bu\\y, + c(£ 7 0 w, £ 7011 ) 

= ||Su|| 2 v + c(£ 7 oit, £ 70 u — u) + c(£ 7 ow — u, u) + c(u, u ) 

< \\Bu\\y, + 11£ 70 "tt - -u||[/|M|t/ + c(u,u) 

<(l + d 1 )\\Bu\\y, + 5\\u\\fj + c(u,u). 

Subtracting the last term for a sufficiently small <5 > 0 this proves the existence of /3 q > 0 such 
that, for all j3 > (3o, there holds ||m||^ < /3\\Bu\\y, + c(u,u) = b(u, Qpu) + c(u,u). This finishes 
the proof of Theorem 0 □ 

3.3 Proof of Theorem [8] 

The proof follows the same lines as the proof of Theorem [ 6 ] in Section 13.21 The only difference 
is that we have to prove ellipticity of c(-, •) = c s i(■, •) on the kernel of B , cf. (LSD . To see this, we 
could again make use of the identity V 7 n £ 7 oU = (^ +/C) 7 o u from Lemma 1121 However, applying 
techniques from |341 135] , our proof extends to more general operators in (fTa|) . 
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Let us abbreviate <f> := 7 n £7o u and let u £ J£ 1 (Q) be the first component of £70 u. Then 
(0,S) r = ||VS||2 2(n)) and we calculate 

c(£'y 0 u,£'y 0 u) = in n £^ou , V^ n £^ou + (| - K)u)r + |(1, V(7£7 0 'u))r| 2 

= (0,V0) r - ((5 + K!)^>,u)t + ( 4 >,u ) r + |(1 , V( 7 £ 7 0 u))r| 2 
= WbV0) r - ((3 + /C')0,«)r + l|Vu||| 2(n) + |(1 ,¥( 7 £ 7 0 w))r| 2 - 

By using that d n Vil> = (| + K.')^ for i/j £ iL _1 / 2 (T), it can be shown that 

-((1 + /C>,£)r > -i^,V^)r - £||VS||£ a(n) . (29) 

For simplicity we refer to [2] Proof of Theorem 9] for a proof. By the same reference, [El, 
Lemma 10], there holds the norm equivalence 

WbWOr + H vs llL 2 (fi) + l( 1 >V(7£7o'it))r| J ^ \\u\\m(n) + U\\ 2 H i/ 2 {r) 

so that 

c{£-f 0 u,£'y 0 u) > i(^,V<^)r + 3IIVw|||, 2(n) + |(1, V(7£7 0 u))r| 2 > IN|^i (n) . 

It follows that c(£ r jQU,£'you) > ||£7o u llir> which finishes the proof. □ 

3.4 Proof of Theorem [9] 

As for the proof of Theorems [6] and |H1 in Sections 13.21 and 13.31 respectively, we only have to show 
the [7-ellipticity of the bilinear form b(-, ®p-) + c ca (-, •). We now make use of Lemma fl3l instead 
of Lemma E2 As in Section T3.21 we have 

\\u\\u < II U - £^ 0 u\\ij + W^ouWf; < \\Bu\\lr, + ||7 0 w||^i /2(r) . 

Applying Lemma 1131 we deduce that 

1170 w ||^fl/2(p) fs ||70 w ||#l/2( r ) + ||7n w ||^-l/2( r ) 

~ (W’7o'u,7oi0r + (7nW,V7 n u)r + |(l,V7u)r| 2 

= (W70W + (3 + /C')7nw,7 0 ii)r + (7nw, (3 - /C)7 0 u + V7 n w)r + |(l,V 7'u)r| 2 - (7n«,7o«)r 
= c ca (u, u) - (7n (u - £70 u ), 7 0 u)r - (7n£7o^, 7o' u )r- 

Since £70 is a harmonic function there holds (7n£7o M , 7o' u )r = (Vu, Vu)q > 0 with u being the 
first component of £70 u £ U. We end up with 

ll7o«||^ 1/2(r) < c ca (u,u ) - (7n(u -£7 0 u) ,7 0 u)r < c ca (u,u ) + ||u - £7 0 u||£/|M|f/- 

To hnish the proof, we argue as in the proof of Theorem El i.e., applying Young’s inequality and 
the continuity ||tt — £^qu\\jj < ||v"'- □ 
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4 Numerical experiments 


In this section we present numerical experiments in two dimensions for the three coupling schemes 
from Sections 12.5.11 (with least-squares BEM), 12.5.21 (with hypersingular operator), and 12.5.31 
(with weakly singular operator). We also discuss their implementation. Numerical results for 
the coupling with both Calderon equations from Section 12.5.41 are very similar to the ones of the 
other variants and are not reported. 

Throughout, we consider regular and quasi-uniform triangulations with compact triangles. 
For uniform refinements we have a sequence of uniform meshes To, ■■■ ,Tl ='■ T of 12. Here, 
Tt is the uniform refinement of 7^_i, i.e., every triangle from 7^_i is split into four triangles. 
Moreover, Si denotes the skeleton induced by Tt and, as before, S is the skeleton of T■ The 
mesh-sizes ho > • • • > =: h denote the largest diameters of the elements in Te (l = 0,..., L). 

Note that there holds hi = 2~^ho- 

Let P p (Ti) be the space of 7^-elementwise polynomials of degree < p, and S 1 (Si) the space of 
edge-wise affine and globally continuous functions on Si . Furthermore, P°(Si) denotes the space 
of edge-wise constant functions on Si. Then we choose the lowest-order discrete trial spaces 

U t := P\Ti ) x [P°(Ti)] d x S\S e ) x P°(Si) 

and our ansatz space for all variants is Uh p := Ul- 

The implementation of the trial-to-test operator 0^ : U — > V requires an approximation. A 
standard way is to consider an enriched subspace V C V that uses the same mesh T as Uhp but 
is enriched by increasing polynomial degrees. We select V := P 2 (T) X [P 2 (T)]~ and consider, 
instead of 0^, the approximated trial-to-test operator 

0: Uhp^-V: (Qw ,v)y = b(w,v) € R. 

Note that we have selected the weight j3 = 1 and, in fact, all our numerical results are based on 
this choice. For an analysis of the approximation of the trial-to-test operator we refer to [ 23] . 
There, it is shown that the increase of polynomial degrees by two in two dimensions (to generate 
the enriched space V) maintains validity of the inf-sup condition. 

The boundary integral operators appearing in the Galerkin bilinear forms of our coupling 
schemes are implemented in a standard way. However, implementation of the least-squares 
bilinear forms involving boundary integral operators is more complicated. This is discussed in the 
next section. Throughout, our examples are conducted in MATLAB and for the discretization 
of boundary integral operators we use the library HILBERT [1|. 

4.1 Implementation of the coupling with least-squares BEM 

An implementation of our least-squares coupling scheme faces two problems: both the operator 
V(y(-)) € H l / 2 {T) and the il 1,/2 (r)-inner product must be approximated. 

First we discuss an approximation of V( 7 (-)). Let n i : .^(T) —»• S' 1 (5^|r) denote the L 2 - 
orthogonal projection. Here, S^iS^r) is the restriction of S 1 ^) to T. We set := and 
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replace the bilinear form 


b(u,Qiw) + c LS {u,w) = (Bu,Bw) v , + (¥(7(11)), V(7(tn)))^i /2(r) 

(cf. ( 11111 . (1141) 1 by the discretized bilinear form 

a h (u, w ) := (Bu , Bw)v + (n^ ¥(7(11)), 11 ^ ¥(7(w)))^i/ 2(r) Vu, w G £ 4 p . 
In the same manner, we replace the right-hand side functional 

L v (@ iw) + L ls (w) = (L v , Bw)v + (¥(u 0 ,</>o),V(7(in))) iP i/2 (r) 
by 

L h {w) := (L v ,Bw) v/ + (II/ l ¥(uo,^o),n/ l ¥(7(to))) jff i/2( r ) Vw G C 4 p . 
Then, our discrete DPG scheme coupled with least-squares boundary elements is 


v-hp £ Uhp . CLh^UhpjW ) — Lh(w) Vin G Uhp. (30) 

As previously mentioned, we also approximate the trial-to-test operator but this is not analyzed 
here. Quasi-optimality of the discrete scheme is maintained, as the following theorem shows. 

Theorem 14. Let Ph : H l / 2 (Y) —> S' 1 (<S|r) denote any bounded projection, i.e., there is a 
T-independent constant C > 0 such that 

\\P h u\\ H i/ 2 {r) < C , ||«|| H i /2(r) Vn G T 1 / 2 (r) and P h u h = u h \/u h e S 1 (S\ r ). 

Then, \ah(u,w)\ < || w ||[/|| ,y \\u f or u,w £U and 

IMIf/ - \\Bu\\ V f + ||P /l ¥( 7 u)|| H i / 2 (r) Vw G U hp . (31) 

The involved constants do not depend on T. 

In particular, one can choose the L^-orthogonal projection Ph = II/,. 

Proof. First we note that the boundedness of Ph implies boundedness of |a^(-, -)|. Furthermore, 
the boundedness of Ph and stability of the continuous problem imply 

\\ Bu \\v + \\Ph ¥( 714 ) 11 ^ 1 / 2 ^) < \\Bu\\v + II V( 7 'u)|| jf/ i/ 2 (r ) ~ IMIc/ Vu G U h p. 

The remainder of the proof is a slight modification of the proof of Theorem 01 By the projection 
property of Ph we have 

\\ u \\u 1$ \\Bu\\v + ||7o M ll// 1 / 2 (r) = \\Bu\\ v > + \\Ph'yo u \\H l / 2 (T) Vu G Uhp- 
Relation 70 Uh p = ¥(7 £ 70 Uh p ) by Lemma [T2] leads to 
||P/i7ow||#i/2(r) = ||-P/,¥(7f 7 0 u)|| H i/2 (r) < ||P/ l ¥(7u)|| H i / 2 (r) + \\P h ¥(7 5 70 ™ - 7tt)|| H i/2 (r) 

< ||P/ l ¥(7u)|| jH -i/2 (r ) + ||u - £^qu\\u Mu g U hp . 

Hence the bound \\u — £^qu\\u < ||Pw||y' hnishes the proof of (1311) . 

Finally, boundedness of the ^-projection in H l / 2 (T) is discussed, for instance, in [29] and 
the references given there. □ 
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Corollary 15. Let the assumptions of Theorem^l^be satisfied. The discrete problem (13011 admits 
a unique solution u^p € Uh p and, with u € U being the solution of (|14[) . there holds 

\\u-u hp \\jj< inf \\u-w\\u. 

■w&Uhp 


Proof. By Theorem 1141 we can apply the Lax-Milgram lemma. It proves that 
solution. The second Strang lemma (see, e.g., H3 Lemma 2.25]) shows that 


has a unique 


1“ - u h P \\u < inf 
weu hp 


I U — w\\u + sup 

■weu hp 


\a h (u,w) - L h (w )| 


\w\\u 


The definitions of •) and Lf, infer that 

a h {u,w) - L h (w) = (Bu — L, Bw) v , + (. P h Vfi/u - (u 0 , fio)), P h V(7™))#i/2(r)- 

Since u is the exact (unique) solution, we know by Theorem [4] that it satisfies equations (fT2l) 
and m, that is, Bu = L and V(y u) = Y(uq, 4>q). □ 

It remains to discuss some implementational aspects. Let {uj- j = 1,..., dim(C7/ ip )} and 
{r/j; j = 1,..., dim(S' 1 (<S^|r))} denote, respectively, bases of Uh p and S ,1 («5>£|r)- Define the 
matrices 


O jk ■= {Vj ,V(7W fc ))r, M jk := (r]j,r]k) r- 

A matrix representation (of the coefficients in the corresponding bases) of II/ l V( 7 (-)) is then 
given by M _1 0. In order to replace the Lf 1 / 2 (T)-inner product for functions Uh G 5' 1 (5^|r) let 
P denote a symmetric and positive definite matrix such that 

x 7 Px ~ ||rf/i||^i/ 2 ( r ) Vrtfe G 5 1 («Sg|r) with coefficient vector x. 

Let x, y denote the coefficient vectors corresponding to Uh,Wh € 5 1 (5£|r), respectively. Then, 
instead of the evaluation of (uh , v h)H 1 / 2 {Y) > we use y T F* x 1 11 particular, we replace the evaluation 
of (TI h V('y(u hp )) ,n h 'y(w hp )) Hl/2 ^ with its matrix form 

y T (0 T M _1 PM _1 0x). 

Different choices of P are possible. For instance, it is well known that the following stabilization 
of the hypersingular integral operator satisfies 

(Wuh,uh)r + | (uh, l)r| 2 — ||«h||#i/ 2 (r) G ^(S^r)- 

We define the corresponding inner product as (uh , Vh)w := (Wu^ , Vh )r + (uh , l)r(^/i j l)r- 

Another possibility is to use multilevel norms. By [33J Theorem 1] and references given there, 
we have the equivalence 

OO 

H w llHi/ 2 (r) ^E 2£ ll^- n ^) u llL(r) Vu € iL 1/2 (T). 

£=0 
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Here. n_i := 0 and is the L 2 -orthogonal projection onto <S ll («S^|r) for l > 0 and, as previously 
defined, 7o,71, ... is a sequence of uniformly refined triangulations. Obviously, if u is discrete, 
the above sum is finite. In particular, we have 

L 

Klltfi/ 2(r) - 2<? IK© - ^-i)uh\\l 2 (r) Vu h £ ^(5 L | r ) 

e=o 

with corresponding multilevel inner product 

L 

(uh,vi ,)ml := ^2 £ ((n^ -Ui_i)u h ,v h )r Vu h ,v h £ S^S^Ir). 
e=o 

Then, (u h ,u h ) M L - IKII^i/ 2 (r) for all u h £ 

Remark 16. (i) Theorem [7^] and Corollary [75] remain valid if we replace the H l / 2 (T)-inner 
product in the definition ofah(-,-) and L/j(-) (•,‘)w or (-,-)ml- 

(ii) The Galerkin matrix ofah(-, •) is symmetric and positive definite. Thus, standard iterative 
schemes such as the CG algorithm can be used. 

(in) Since the bilinear forms Ch y ,c s i,c ca are not symmetric, it follows that the Galerkin ma¬ 
trices of fe(-,0 j g(-)) + c(-, •) with c £ {ch y , c s i, c ca } are non-symmetric. However, their symmetric 
parts are positive definite and the GMRES method can be used as solver. 


4.2 A simple error estimation 

In our implementation, we approximate the data rto,^o by their ^-projections uqi l € 5 1 (<S|r), 
(J)q h £ R°(5|r)- Then we can use the discretized boundary integral operators for the evaluation 
of the right-hand sides. Clearly, this induces additional consistency errors (usually called oscil¬ 
lations). They are of standard type and are therefore neglected in the following considerations. 

The next result provides us with a reliable bound for the distance of any discrete function to 
the exact solution. 


Theorem 17. Let the assumptions of Theorem If be satisfied, let u^h £ («S|r), <ftoh £ R°(5|r) 
be approximations of the given data uo,(f>o, and let u be the exact solution of problem © as 
defined in Theorem^ with data uq,4>o being replaced by uq^, <poh- Then, 

\\u-w hp \\u < \\Bw hp - L\\ v > + \\Ph^{(u 0 h,(foh) -~T w hp)\\ H i/2(r) 

+ \\h 1/2 \7 r H{(u 0h ,<j) 0h ) -7«bi P )||L 2 (r) £ U hp 

with T-independent constants. 

Proof. The boundedness and t/-ellipticity of 6(-, ©i •) + cls(' 5 •) (cf. (12711 ) show 


u - w hp \\u ~ || B(w hp - u)\\ V t + ||¥(7U - 7W/ip)IIh 1 /2(d 

< \\B(w hp - u)\\ v > + WPnYfiyu - 'yw hp )\\ H i/2 ( r ) 
+ 11(1 - V(7u -7ii>/, p )|| H i/2 (r) . 
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Using the fact that u is the exact solution, we have Bu = L and V( 7 u) = H(uoh,4>oh)- Then, 
applying |18( Lemma 24], that is, 

||(1 - P h )^{(u oh ,(p oh ) - -/Whp)\\ H i/^r) ^ ||h 1/2 VrV((u o/l ,0o h ) - 7™fcp)ll£, 2 (r). 

we finish the proof. □ 

From the last theorem we infer that, up to data oscillation terms, we have a computable error 
bound 

est(u^p) . est(ii/ 1 p, f, UQh, (f>oh) • \\Buh p L||y/ T 1111/^ V((rto/i, (froh) 'T^ , hp)\\H 1 / 2 (r) 

+ ||h 1//2 Vr y((u 0 h,(j)oh) ~ l u h P )\\ L2 (T)- (32) 

In particular, the first term is the usual energy error in the DPG method and the second term can 
be evaluated as described in Section IQ above. The last term has the same form as a weighted 
residual error estimator in BEM, see |l8j for an overview. Note that est(-) can be used to estimate 
both the error of the least-squares coupling and of the coupling methods with Galerkin BEM. 

Remark 18. For the coupling with Galerkin boundary elements (at least for the case with hy¬ 
persingular operator from Subsection \2. 5. 2) . reliable and localizable residual error estimators can 
be developed, following the lines ofW and references given there. However, for brevity and better 
comparability of the least-squares and Galerkin boundary element couplings, we use est(-) from 
(13211 to plot an upper bound for the global error (comprising all solution components) for all our 
coupling methods. 


4.3 Example with smooth solution 

For the first example we take the L-shaped domain U sketched in Figure Q] and consider the 
transmission problem (JT]) with prescribed smooth solution 

u(x, y) = - , u c (x, y ) = 0. 


The corresponding data are / = —2, uq = u |r, 4>0 = <9 n n u. The definition of the trace norms 
II ‘ |li : f i / 2 ( 5 ), || ■ ||ff-i/ 2 ( 5 )j the quasi-optimality of the methods and standard approximation theory, 
cf. [5] US] , yield the a priori estimate 


\\u-u hp \\u< inf \\u-w hp \\u < /i s (|M| H i+ s(n) + H/lltf^n)), s<l. 

Hence, we expect a convergence order 0[h). This order is confirmed by Figure Q] for the least- 
squares method (right plot), and for both one-equation coupling methods with Galerkin boundary 
elements (Figure [2]). Here, for the least-squares BEM part, we have used the multilevel inner 
product (•, - )ml- 
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Figure 1: Domain with initial triangulation 7o (left); L 2 errors of u and er, and error estimator 
est(-) for the least-squares coupling (right), with smooth solution u. 



Figure 2: L 2 errors of u and a, and estimator est(-) for coupling with Galerkin BEM and smooth 
solution u: c(-, •) = Ch y (-, ■) (left) and c(-, ■) = c s i(-, •) (right). 
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Figure 3: L 2 errors of u and <x, and estimator est(-) for least-squares coupling and singular 
solution: uniform meshes (left) and adaptively refined meshes (right). 




Figure 4: L 2 errors of u and er, and estimator est(-) for coupling with Galerkin boundary elements 
and c(-, ■) = Ch y (-, •) with singular solution: uniform meshes (left) and adaptively refined meshes 
(right). 


23 







Figure 5: L 2 errors of u and er, and estimator est(-) for coupling with Galerkin boundary elements 
and c(-, •) = Csi(-, •) with singular solution: uniform meshes (left) and adaptively refined meshes 
(right). 

4.4 Example with singular solution 

We consider the L-shaped domain from Figure |T| and prescribe the exact solution 

u(x, y) = r 2/3 cos(2/30), u c (x , y) = ■ - + ^, 8 2 ■ 

10 [x-^Y + y 2 

Here, (r, 9) are polar coordinates centered at the incoming corner. There holds u G H S (Q) for 
all s < 5/3 and, therefore, we expect an overall convergence rate of 0{h , 2 / 3 ) = 0(N ~ 1 / 3 ). This 
is confirmed for all coupling methods by the left plots of Figures [31 HI El We observe that the L 2 
errors of u converge with a better rate between 0(N~ 5 / 12 ) and 0(N ^ 1 / 3 ). 

For this example we used (•,-)w as inner product in T). This is the simplest selec¬ 

tion that also works for non-uniform meshes, needed when considering an adaptive refinement 
strategy. For such a strategy we use the quantities 

II L - Bu\\y, + H/1 1 / 2 V r V((uo/i, 4>oh) ~ 7 «fc P )lli 2 ( r) 

as indicators, and refine with bulk-strategy (also called Doerfler-marking) and 30 percent marking 
(of the squared quantities). This error estimator is purely heuristical and is used for its locality. 
From Figures El HI and E] (right plots) we observe that the rate 0(N ~ 1 / 2 ) is recovered for all 
coupling methods. 
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